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Abstract 

Lorentz-violating operators involving Standard Model fields are tightly constrained by exper- 
imental data. However, bounds are more model-independent for Lorentz violation appearing in 
purely gravitational couplings. The spontaneous breaking of Lorentz invariance by the vacuum 
expectation value of a vector field selects a universal rest frame. This affects the propagation of 
the graviton, leading to a modification of Newton's law of gravity. We compute the size of the 
long-range preferred-frame effect in terms of the coefficients of the two-derivative operators in the 
low-energy effective theory that involves only the graviton and the Goldstone bosons. 
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The possibility of breaking Lorentz invariance is of interest to theoretical physicists for a 
variety of reasons. For instance, spontaneous Lorentz violation has been proposed at various 
' ;imes as a starting point in alternative theories of electrodynamics and of linear gravity 
3], and as a possible solution to the horizon problem in cosmology [^. Lorentz violation 
has also been discussed as a potentially observable signal of physics beyond the Planck 
scale (whether in the context of string theory, noncommutative geometry, or loop quantum 
gravity) 0], and some researchers have claimed that there is evidence of such violation in 
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the measurements of the energies of cosmic rays 

Experimental data puts very tight constraints on Lorentz violating operators that involve 
Standard Model particles lal, but the bounds are more model-independent on Lorentz vio- 
lation that appears only in couplings to gravity |7|,|8|]. One broad class of Lorentz-breaking 
gravitational theories are the so-called vector-tensor theories in which the space-time metric 
g^ v is coupled to a vector field which does not vanish in the vacuum. Consideration of 
such theories dates back to [2| and their potentially observable consequences are extensively 
discussed in . These theories have an unconstrained vector- field coupled to gravity. The- 
ories with a unit constraint on the vector field were proposed as a means of alleviating the 
difficulties that plagued the original unconstrained theories . 

The phenomenology of these theories with the unit constraint has has been recently 
explored. It has been proposed as a toy model for modifying dispersion relations at high 
energy The spectrum of long- wavelength excitations is discussed in jl^, where it was 
found that all polarizations have a relativistic dispersion relation, but travel with different 
velocities. Applications of these theories to cosmology have been considered in ^ 
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Constraints on these theories are weak, as for instance, there are no corrections to the Post- 
Newtonian parameters 7 and (5 Jig]. The status of this class of theories, also known as 
'Aether-theories', is reviewed in |l7| . 

Here we begin by considering the general low-energy effective action for a theory in 
which Lorentz invariance is spontaneously broken by the vacuum expectation value (vev) of 
a Lorentz four- vector S^. With an appropriate rescaling, the vev satisfies 

(S»S») = 1 , (1) 

since we assume the vev of is time-like. The existence of this vev implies that there 
exists a universal rest frame (which we sometimes refer to as the 'preferred frame') in which 
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S 11 = 5q . When the resulting low-energy effective action is minimally coupled to gravity, we 
shall see that it simply becomes the vector-tensor theory with the unit constraint. 

Objects of mass Mi and M 2 in a system moving relative to the preferred-frame can 
experience a modification to Newton's law of gravity of the form 

HQ 

^Newton - ~<^N " I 1 ~ ~ ~ 2 I {2.) 

where w is the velocity of the system under consideration, such as the solar-system or 
Milky Way galaxy, relative to the universal rest frame. The main purpose of this note is 
to compute a 2 in theories where Lorentz invariance is spontaneously broken by the vev of 
a four-vector. This PPN coefficient is more strongly constrained by experiment than the 
other PPN parameters 7 and (3 [3], so it is natural to focus on it. 

The vev of spontaneously breaks Lorentz invariance. But as rotational invariance is 
preserved in the preferred frame, only the three boost generators of the Lorentz symmetry 
are spontaneously broken. The low-energy fluctuations S^^x) which preserve Eq. ((TJ are 
the Goldstone bosons of this breaking, i.e., those that satisfy 

S^S^x) = 1 . (3) 

In the preferred-frame the fluctuations can be parameterized as a local Lorentz transforma- 
tion 

S*{z) = Aft(x) = -r=L= ( V| . (4) 
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Under Lorentz transformations S^^x) — > h.^S u (x) and the symmetry is realized non- 
linearly on the fields <p l . Using this field S ,J '(x) we may then couple the Goldstone bosons 
to Standard Model fields. Since however, the constraints on Lorentz-violating operators 1 
involving Standard Model fields are considerable |6jij, we instead focus on their couplings 
to gravity, which are more model independent because they are always present once the 
Goldstone bosons are made dynamical. 

The Goldstone bosons are made dynamical by adding in kinetic terms for them. Since 
Lorentz invariance is only broken spontaneously, the action for the kinetic terms should 



1 More correctly, operators that appear to be Lorentz violating when the Goldstone bosons <j? are set to 
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still be invariant under Lorentz transformations. The only interactions relevant at the two 
derivative-level and not eliminated by the constraint Eq. Q are 

£ = c 1 d a S f3 d a Sp + (ca + c z )d^d v S u + c i S>*d pj S a S v d 1/ S a . (5) 

Expanding this action to quadratic order in 0*, one finds that the four parameters Cj can be 
chosen to avoid the appearance of any ghosts. In particular, we require c\ + c 4 < 0. 

With gravity present the situation is more subtle. One expects t he g ravitons to 'eat' 
the Goldstone bosons, producing a more complicated spectrum [3, The covariant 

generalization of the constraint equation becomes 



g^(x)S^x)S u (x) = 1 (6) 

and in the action for we replace d p — > V M . Note that there is no "Higgs mechanism" to 
give the graviton a mass, since the connection is linear in derivatives of the metric. 

Local diffeomorphisms can now be used to gauge away the Goldstone bosons. For under 
a local diffeomorphism (which preserves the constraint Eq. (JHJ)), 

S*V) = ^{x) (7) 

and with x'^ = + e M , = v M + 

^(x 1 ) = <jf{x) + v p d p e^ (8) 

from which we can determine e M to completely remove M . Note that in the preferred frame, 
e l can be used to remove 0\ In this gauge, the constraint Eq. @ reduces to 

S°(x) = (l-h 00 (x)/2) . (9) 

The residual gauge invariance left in e° can be used to remove hoo. This is an inconvenient 
choice when the sources are static. In a more general frame with (S^) = v M , obtained by a 
uniform Lorentz-boost from the preferred-frame, the constraint Eq. (jUJ) is solved by 

S»(x) = (1 - v p v a h pa {x)/2) . (10) 

Next we discuss a toy model that provides an example of a more complete theory, that 
at low-energies reduces to the theory described above with the vector field satisfying a unit 
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covariant constraint (jH}- 2 Consider the following non- gauge invariant theory for a vector 
boson A 11 , 

C = -\g» v g pa V p A»V a A v + A (g^A»A v - v 2 f . (11) 
Fluctuations about the minimum are given by 

g^u = vv» + K» > A>1 = vfl + V ■ ( 12 ) 

This theory has one massive state $ with mass M$ oc A 1//2 t> , which is 

$ = + . (13) 

In the limit that A — > oo this state decouples from the remaining massless states. In the 
preferred frame the only massless states are h^, and Since we have decoupled the heavy 
state, we should expand 

A = v + \f + vh m /2\ - vh 00 /2 -+v- vh 00 /2 (14) 

where in the last limit we have decoupled the heavy state. Note that this parameterization 
of A is precisely the same parameterization that we had above for S°. In other words, in the 
limit that we decouple the only heavy state in this model, the field A^ satisfies g^ v A^A u = v 2 , 
which is the same as the constraint (JHJ) with A^ — > vS^. 

In the unitary gauge with <p l = 0, the only massless degrees of freedom are the gravitons. 
There are the two helicity modes which in the Lorentz-invariant limit correspond to the two 
spin-2 gravitons, along with three more helicities that are the Goldstone bosons, for a total 
of five. The sixth would-be helicity mode is gauged away by the remaining residual gauge 
invariance. 

But the model that we started from does have a ghost, since we wrote a kinetic term 
for N 1 that does not correspond to the conventional Maxwell kinetic action. The ghost in 
the theory is A , which in our case is massive. The presence of this ghost means that this 
field theory model is not a good high-energy completion for the low-energy theory involving 
only and gravity which we are considering in this letter. We assume that a sensible high 
energy completion exists for generic values of the q's. 



For a related example, see [ipf . 
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Now we proceed to compute the preferred-frame coefficient «2 appearing in the modifi- 
cation to Newton's law. 
The action we consider is 

S = Jd A x^ (£ EH + C y + £gf) (15) 

with 3 

£ EH = -^R (16) 

C v = dVaS^Sp + c 2 W^W v S v + cgV^V^ + c 4 ^V „S a S v V v S a , (17) 

and we use the metric signature (H ). This is the most general action involving two 

derivatives acting on S** that contributes to the two-point function. Note that a coefficient 
C3 appears, since in curved spacetime covariant derivatives do not commute. Other terms 
involving two derivatives acting on may be added to the action, but they are either 
equivalent to a combination of the operators already present (such as adding R^S^S"), or 
they vanish because of the constraint Eq. We assume generic values for the coefficients 
Cj that in the low energy effective theory give no ghosts or gradient instabilities. 

As previously discussed, satisfies the constraint (0). We also assume that it does not 
directly couple to Standard Model fields. In the literature, Eq. © is enforced by introducing 
a Lagrange-multiplier into the action. Here we enforce the constraint by directly solving for 
S^, as given by Eq. (fTUj) . and then insert that solution back into the action to obtain an 
effective action for the metric. 

In our approach there is a residual gauge invariance which in the preferred-frame corre- 
sponds to reparameterizations involving e° only. To completely fix the gauge we add the 
gauge-fixing term 

£ gf = -| (S p S a S' 1 d^h pa ) 2 . (18) 
Neglecting interaction terms, in the preferred frame the gauge-fixing term reduces to 

£ g f = ~ (doh 00 ) 2 ■ (19) 



3 The coefficients Ci appearing here are related to those appearing in, for example 1 1 3j . by c^ crc 

_ c there/ 167rG f 
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Physically, this corresponds in the a — > oo limit to removing all time-dependence in hoo, 
without removing the static part which is the gravitational potential. This is a convenient 
gauge in which to compute when the sources are static. 

At the two- derivative level, the only effect in this gauge of the new operators is to modify 
the kinetic terms for the graviton. The dispersion relation for the five helicities will be of 
the relativistic form E = 8\k\, but where the velocities 8 are not the same for all helicities 

n 

and depend on the parameters Cj |13||. This spectrum is different than that which is found 
in the 'ghost condensate' theory, where in addition to the two massless graviton helicities, 
here exists a massless scalar degree-of-freedom with a non-relativistic dispersion relation 
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There exists a range for the q's in which the theory has no ghosts and no gradient 
instabilities Q]. In particular, for small q's, no gradient instabilities appear if 

Ci + c 2 + c 3 Ci 

> and > . (20) 

Ci + C 4 Ci + c 4 

The condition for having no ghosts is simply c\ + c 4 < 0. 

The correction to Newton's law in Eq. (j2J) is linear order in the source. Thus to determine 
its size we only need to find the graviton propagator, since the non-linearity of gravity 
contributes at higher order in the source. In order to compute that term we have to specify 
a coordinate system, of which there are two natural choices. In the universal rest frame the 
sources, such as the solar system or Milky Way galaxy, will be moving and the computation 
is involved. We instead choose to compute in the rest frame of the source, which is moving at 
a speed \w\ -C 1 relative to the universal rest frame. Observers in that frame will observe the 
Lorentz breaking vev ~ (1, —w). In the rest frame of the source, a modified gravitational 
potential will be generated. Technically this is because terms in the graviton propagator 
v ■ k ~ w ■ k are non-vanishing. It is natural to assume that dynamical effects align the 
universal rest frame where = Sq with the rest frame of the cosmic microwave background. 

In a general coordinate system moving at a constant speed with respect to the universal 
frame the Lorentz-breaking vev will be a general time-like vector v M . Thus we need to deter- 
mine the graviton propagator for a general time- like constant v M . Since Lorentz invariance is 
spontaneously broken, the numerator of the graviton propagator is the most general tensor 
constructed out of the vectors v^, k v and the tensor rj pcr . There are 14 such tensors. Writing 
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the action for the gravitons as 



S = ±J d 4 kh a ?(-k)K a p Wp (k)h° p (k) 



(21) 



it is a straightforward exercise to determine the graviton propagator V by solving 

K a ^(k)V^(k) = \ (rj p tf + rfoft . (22) 

The above set of conditions leads to 21 linear equations which determine the 14 coefficients 
of the graviton propagator in terms of the coefficients q and the vev v^. Seven equations 
are redundant and provide a non-trivial consistency check on our calculation. 

Although it is necessary to compute all 14 coefficients in order to invert the propagator, 
here we present only those which modify Newton's law as described previously (assuming 
stress-tensors are conserved for sources). These are 

P£ = {Arf V CT + B(?f V 7 + ?f V P ) + C(v a v fi rT + v p v a r] ap ) 

+Dv a v l3 v p v a + E(//^ + v a v a rf p + v^v p rf c + v (3 v a r) ap ) } (23) 

We find that each of these coefficients is independent of the gauge parameter a. We also 
numerically checked that without the presence of the gauge-fixing term the propagator could 
not be inverted. 

To compute the preferred- frame effect coefficient a 2 , we only need to focus on terms in 
the momentum-space propagator proportional to (v ■ k) 2 . To leading non-trivial order in 
G(v ■ k) 2 and in the q's we obtain, from the linear combination A + IB + 2C + D + 4E, 

1 



g 00 = 1 + 8nG 



2cf + Ac 2 3 (c 2 + c 3 ) + 



(27r) 4 P[ 1N k 2 c 1 (c 1 + c 2 + C3 ) 

+c?(3c 2 + 5c 3 + 3c 4 ) + ci((6c 3 - c 4 )(c 3 + c 4 ) + c 2 (6c 3 + c 4 ))]| f 00 (k) (24) 

where in the first line k is a four- vector. Next we use v p = (1, —w), place the source at the 
origin, substitute T 00 = M<5 (3) (£) or T 00 (k) = 2irM5(k°) and use 

1 



r cn 
J (27 



(i 3 /c kjkj g. s _ 

—tC- 



to obtain 



M 



9oo 



1-2G N — 1 



(2tt)3 fc 4 
(w • r) 2 



8nr 



r XiXj 

5 V - 



(25) 



8ttG> 



r 2 2ci(ci + c 2 + c 3 ) 



2c 3 + 4c 2 (c 2 + c 3 ) + 



+c 2 (3c 2 + 5c 3 + 3c 4 ) + ci((6c 3 - c 4 )(c 3 + c 4 ) + c 2 (6c 3 + c 4 )) 



(26) 



where we have only written those terms that give a correction to Newton's law proportional 
to [w ■ f/r] 2 . We have also assumed that \w\ C 1 so that higher powers in w ■ r/r can 
be neglected. The factor of 1/ci in the preferred-frame correction to the metric arises 
because when c\ —>■ the "transverse" components of 0* have no spatial gradient kinetic 
term. Similarly, the factor of l/(ci + c 2 + c 3 ) arises because when c\ + c 2 + c 3 — > the 
"longitudinal" component of tp 1 has no spatial gradient kinetic term. Either of these cases 
causes a divergence in the static limit. 4 

The coefficients q redefine Newton's constant measured in solar system experiments and 
we find that 

Gk = G [1^8 tG (c 1 + c 4 )1^ 1+8 ^ Ci+C4) (27) 

which agrees with previous computations to linear order in the Cj's after correcting for the 
differences in notation [l^. Il7|. 

The experimental bounds on deviations from Einstein gravity in the presence of a source 
are usually expressed as constraints on the metric perturbation. Since the metric is not 
gauge-invariant, these bounds are meaningful only once a gauge is specified. In the litera- 
ture, the bounds are typically quoted in harmonic gauge. Here, the preferred-frame effect is 
a particular term appearing in the solution for Hqq. For static sources, the gauge transfor- 
mation needed to translate the solution in our gauge to the harmonic gauge is itself static. 
But since a static gauge transformation cannot change /ioo ; we may read off the coefficient 
of the preferred-frame effect in the gauge that we used. 

By inspection 

« 2 = 7 8 ^ Gn , r [2c? + 4c|(c 2 + c 3 ) + c?(3c 2 + 5c 3 + 3c 4 ) 

Cl(Ci + c 2 + c 3 ) 

+c 1 ((6c 3 -c 4 )(c 3 + C4) + c 2 (6c 3 + c 4 ))] , (28) 

which can be compared with the experimental bound |a 2 | < 4 x 10~ 7 given in [18]. 

A considerably stronger constraint on the size of the q's can probably be obtained from 
the gravitational Cherenkov radiation of the highest-energy cosmic rays 
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4 This divergence can of course be avoided by considering higher-derivative terms in the action for the 
Goldstone bosons. This would then give non-relativistic dispersion relations for these modes (as was the 
case in |2l|). 
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